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Abstract. - We point out the existence of finite charge and spin Hall conductivities of graphene 
in the presence of a spin orbit interaction (SOI) and localized magnetic impurities. The SOI in 
graphene results in different transverse forces on the two spin channels yielding the spin Hall current. 
The magnetic scatterers act as spin-dependent barriers, and in combination with the SOI effect lead 
to a charge imbalance at the boundaries. As indicated here, the charge and spin Hall effects should 
be observable in graphene by changing the chemical potential close to the gap. 



Introduction. — Graphene, currently under intense experimental and theoretical in- 
vestigations [1-4], consists of a single layer of carbon atoms arranged in a honeycomb lattice 
made of two interpenetrating, cr-bonded triangular sublattices, A and B. The low-energy 
spectrum and the transport properties near the neutrality point are dictated by the 7r and 
7r* bands that form conical valleys touching at the high symmetry points K and K' of the 
Brillouin zone [5]. A key point is the linear dispersion near K and K 1 which renders a low- 
energy theory in terms of an effective Dirac-Weyl Hamiltonian [4] . This results in a number 
of fascinating phenomena such as the half-quantized Hall effect [1,2,6,7]. Remarkably, the 
spin relaxation length is measured to be as long as ~ 1.5 fj,m in low- mobility devices [8] 
which makes graphene interesting for spintronic applications. In this regard the spin-orbit 
coupling (SOI) [9] is a decisive factor as it determines the spin-decoherence time and opens 
the way for the spin manipulation and control. Generally, in ferromagnets SOI plays a key 
role in the anomalous Hall effect (AHE) [10-15] and it is also essential for the spin Hall 
effect (SHE) [16-19]. Manifestations of SOI in graphene have been addressed by Kane et 
al. [20,21] indicating that the spin Hall conductivity in the undoped graphene is quantized 
due to a gap produced by SOIs in the absence of a magnetic field. Following Kane et al.'s 
work, Sinitsyn et al. [22] predicted a substantial spin Hall conductivity in doped graphene 
due to the skew scattering, and find that the charge Hall conductivity vanishes because of 
a cancelation of contributions of bands with opposite spins. 

In this work, we study the Hall effect in graphene in the presence of SOI and localized 
magnetic impurities. We find the charge Hall effect is generally finite as a result of the 
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combined influence of SOI and spin flip scattering at the magnetic impurities, a mechanism 
different from the skew scattering [11] and the side jump [12] mechanisms established for 
AHE in ferromagnetic metals. The SOI in graphene exerts an asymmetric transverse force 
on electrons with different spins, generating thus a spin Hall current. In addition, the 
scattering off the magnetic impurities serves as a spin dependent barrier leading to a charge 
imbalance at the boundary of the sample. We argue that the charge and the spin Hall 
effects should be observable in the presence of SOI and magnetic impurities by changing the 
chemical potential close to the gap. 

Theoretical model. — The Hamiltonian for the clean undoped graphene with SOI 
is [20,23] 

H = ^V k (^k ■ (JT Z + Aa z s z )ip k , (1) 

k 

where 1p k = (c WAK p C kBK ^, C kBK ,^, C kAK ,^, C kAK ^, CkBKl' C kBK'l^ C kAK'l)- c licrTs 1S tne 

creation operator of a single-particle state labeled by the momentum k and associated with 
the sublattice a = A, B and the Dirac point r = K,K', and having the spin s =t, J., vf 
is the Fermi velocity and A is the SOI strength parameter. Upon doping with localized 
magnetic impurities the graphene carriers couple to the impurities via 

^^E^ +q ^(q)^k, (2) 

k,q 

where V is the volume of the system, and 

Vi(q) = ( To ^ U ^P^ + ^ Vs. (3) 

Pi (q) stand for the Fourier components of the impurity density, Q is a vector pointing from 
the point K to the point K', u\ and u 2 are the intra-valley and inter-valley scattering 
strengths, respectively, and S is the spin operator of the localized magnetic impurities [24]. 
The case of a single localized impurity has been treated in [25] . To obtain the spin-dependent 
lifetimes of the quasi-particles and their energetic positions in the presence of impurities we 
calculate the spin-dependent imaginary and real parts of the self-energies E, which in the 
Born approximation [26] reads 

E(k, k'; lo + i V ) = £ - q)S (q; w + ir?)^(q - k')) U) 

q 

where (...) stands for the impurity average and C?o is the resolvent of H. Straightforward 
calculations show that in spin space E has the structure 

The matrix expressions for the impurity averaged, retarded (R) or advanced (A) Green 
functions are [27] 



G fl / A (k,( 



{u-Y,l H/A )?-vl\?-{&+Y,l H/A )-> U 

q u-^ R/A +v F k-aT z -(A-^ A )a z 



\ (u 1 -Y, l a/A r-vlk.i-(&~Y, l c R/A Y J 

(6) 

From iJSJ) it is obvious that E consists of a spin-dependent but sub-lattice independent part 
E , whereas E c depends on both the sublattice index and on the spin. Thereby we find the 
explicit expressions 
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where 



and 



and 



Here, 



r s = 



n,( M 2 + M 2 )(a s + (5 2 )K y-A 2 - D\ 
Ao = ufi~~~A2 ' 

_ 0(0 < lo 2 — A 2 < J 2 )7m,K + u|)(o fl + (5 2 ))M 

_ (-l)^( M 2 -u 2 )(q s -(5 2 ))A r 2 -A 2 - ^ 
47TV 2 1 ^ 2 -A 2 

(-l)'fl(O < uj 2 - A 2 < £> 2 )sgn(o;)7rn i (M 2 - u|)(a a - (5 2 ))A 



a s=T = (5-5+), a s =x = (5+5-), 

(— f) s = ±1 for s =T| 5 an d 5+ = 5a: ± i5 y are the lowering and raising operators of the 
impurity spin. 9{x) is the step function, and rij is the average impurity density. D = vpk c 
and k c is a cutoff for the k summation. 

Charge and spin currents. — The i component of the charge current operator is 

dH 

J % = e—— = ei'FViTz- 
dki 

The spin current with the spin being along the z axis and flowing in the x direction we 
evaluate from the anti-commutator between the velocity operator and the Pauli matrix [19] 

T2 h dH h . 
J * = 4 {s *>dk- } =2 VF<TxTzSz - (?) 

The electric conductivity cry is the (Kubo) linear response function to the external electric 
field [28] 

ImllMi + M;) 
ay = hm , 

where Hf-{u>) is the retarded current-current correlation function obtained from an analytic 
continuation of the Matsubara function 

1 ff 3 

n tf (if2 n ) = -/ dre^iTrJi^JjiO)), fl n = 2wnT 
V Jo 

and (3 is the inverse temperature (T) and T T is the r ordering operator. We calculated cry 
by adopting the ladder approximation for the current vertex determined by Bethe-Salpeter 
equation [29,30]. 

T xx ' (uj, </) = a x r z + ^{Vi(k- k')G A (k', uj)T xx ' (uj, uj')G x ' (k', w ')^(k' - k)>, (8) 

k' 

where A, A' = R, A; uj' — uj + Q and (• ■ ■ ) stand for the average over the impurity distribu- 
tions. By iteration we express Eq.© as (r AA ' = T xx ' (uj, uj')) 

T \\> = a x\' axTz + b XX' ayTz+c xy axTzSz + d XX' ayTzSz ( 9 ) 
Substituting Eq.© in Eq.© we find 

a xx '+c xx ' = l + ao(5_5+)n 2 [S 3 (a AA ' - c xx ' ) + i~ 4 (b xx ' - d xx ' )} + a (S 2 z ) 

n^ia^' + c xx ' ) - iZ 2 (b xx ' + d xx ' )] , (10) 
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a xx ' - c xx ' = 1 + a (S 2 )Il 2 [~ 3 (a xx ' - c xx ') + i~ 4 (b xx ' - d xx ' )] + a {S+S-) 
E x [Si (a x v + c AA ' ) - iS 2 (& AV + d AV )] , 

b xx '+d xx ' = p (S-S+)U 2 [E 3 (b xx ' -d xx ')-i~ 4 (a xx ' -c xx ')]+p (S 2 z ) 
Ei [Hi (fe AA ' + d AA ' ) + zS 2 (a AV + c AA ' )] , 

fe AA'_ rf AA' = ^( J S«)n 2 |S8(6^'-d AA ')-* S 4(o AA '-C AX/ )]+A(^+'S'-> 

Ei [Si (6 AA ' + d AA ' ) + iZ 2 (a AA ' + c AA ' )] , 

wtioro rv — "i("l-"D 13 _ "i("l+"2) TT — TT AA ' " — ~ AA ' anH 

wiiere a - — — , p - — — , i-H = llj , = , and 

(u,-sTU),a ) 2_ (a±s t(;),>. ) 2 

ttAA' _ V-=J a) ' V ) a -(A± = I UM V 

All /n\ 



K 2 ) ( w '-S;< l) ' A ') 2 -(^-S t 1 J <i) ' A )2-(A±S c Ta »' A ')2 + (A±E: a) ' A )2' 

3$, = ( W - sJ (i) ' A )(c' - S^^') — (A ± «(A ± St(i). A '), 
3$, = (A ± £TU), A )(a/ - E^ A ') _ (w _ e tu).a )(a ± E Ta).v }j 

Using Eq. ([9]) we obtain for the Hall conductivity at zero temperature the first central result 

<Jyx = 4? { T i(b AR + d AR ) + T 2 (b AR - d AR ) + T 3 (a AR + c AR ) + T 4 (a AR - c AR )} , (11) 

where a AR , b AR , c AR , d AR are the solutions of Eqs.(fIU |) -([TT1 ) in the limit fl — ► 0, and (/i is 
the chemical potential) 

_ 1 ( A1 -AT(»)^ ( A±AT(i))^ + rT'" 2 -rni) 2 ( A1 -A;»»)^r ( l"» 2 -(A±AT(»)^+rT(»^ 

Jl < 2 ) _ 2 (/1 -Aj">)rJ">±(A ± Al">)rIU> aiCtan 2( M -A^>)r;">±2(A±Al">)rP> ' 

(12) 

_ Atr T(i) ±Ar T(-l) + (A T(-l) r T(^)_ A 1 | (i) r T(l) ) (M _AT^')^-rT'^ 2 -(A±Af^>) 2 +rT(» 2 

J3 ( 4 ) ~ (p-Aju^rJ^'ifAiAl"')^^' arctan 2( M -Aj"))rT(i»±2(A±Al( i ))rI") • 

(13) 

From this procedure we obtain for the spin conductivity 

< x = j^{Ti(a AR + c AR ) - T 2 (a AR - c AR ) - T 3 {b AR + d AR ) + T 4 {b AR - d AR ) + C}, 

(14) 

where C = (1+ao a ° s ( i»tf j ( if f + )(s + S-) ■ The s P in Hal1 conductivity is 

°U = J^i T i( bAR + O - T ^ bAR - dAR ) + T 3(a AR + c AR ) - T 4 {a AR - c AR )}. (15) 

Discussions. — From the above derivation the following general conclusions are in- 
ferred: i) When the SOI strength A = and the impurity spin ladder operators S + , S~ / 0, 
we find for the sublattice-dependent part of the self-energy T,l R ^ A = Y,i R ^ A = 0, the param- 
eters in Ed.([I5]) T 3 = T 4 = 0, and the parameters in Eq.flTUI) S AR = E AR = 0. Substituting 
these results in Eqs.^TDJ-dTTJ, we obtain b AR - d AR = b AR + d AR = 0. This means that the 
charge Hall conductivity a yx = vanishes (see Eq. (fTT)0 as well as the spin Hall conductivity 
a yx ~ ( see Eq-lfTB])). Hence, we conclude that merely spin-flips at the magnetic impurities 
do not induce charge and/or spin Hall effect in absence of SOI. ii) For the SOI strength 
A ^ and the impurity spin operators S + = S~ = 0, our model for the magnetic impu- 
rities is not reducible to the nonmagnetic impurity case; the former case is related to the 
Pauli spin matrix s z , while the latter is expressed in a unit matrix, but they show similar 
behaviour. Because in our case a AR - c AR = a AR + c AR and b AR - d AR = -(b AR + d AR ) 
when the self-energies satisfy the relations: Y,q R ^ A = Sq and Y^ C R ^ A = -So , we 
conclude that the charge Hall conductivity a yx = and the spin Hall conductivity is not 
vanishing; a result which is consistent with the findings of Ref. [22]. The charge Hall con- 
ductivity vanishes because of a cancelation between bands with opposite spins. Hi) For the 
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SOI strength A ^ and the impurity spin operators S + , S~ ^ 0, the charge and spin Hall 
conductivities are in general nonvanishing, similarly as in the case of magnetic impurities 
in two dimensional electron gas with Rashba SOI [30,31]. However, in the latter case, the 
spin Hall conductivity is suppressed to zero by scattering from nonmagnetic impurities with 
a linear Rashba SOI [32]. 

This analysis suggests that the charge Hall effect is the consequence of the combined 
influence of the SOI and the spin flip scattering at the magnetic impurities. This is an es- 
sentially different mechanism from the skew scattering [11] and the side jump mechanism [12] 
in the ferromagnetic metals. In fact, the latter do not occur in our graphene model since 
the SOI is intrinsic and homogenous [33]. Evidence is provided by investigating the longi- 
tudinal spin conductivity. For example, when A = and S + , S~ =/= it is straightforward 
to verify that the spin Hall conductivity vanishes, but a s xx ^ 0, meaning that the scattering 
of spin-up and spin-down electrons by magnetic impurities in graphene is different in the 
longitudinal direction only, but does not yield an effective force in the transverse direction, 
thereby resulting in a non-vanishing spin current and a s yx = 0. In contrast, for A^O and 
S + , S~ = 0, we obtain a* ^ and a s xx — 0. This is due to opposite SOI forces in transverse 
direction on the charge current of the spin-up and spin-down electrons. Thus, we conclude 
that the presence of both the SOI and the magnetic impurity brings about the occurrence 
of the charge Hall effect without external magnetic field in graphene. 

In the Boltzmann limit [34], which in our model is achieved when Ag,A* = and 
H » |r- /c | we find 

4e 2 ((5 + 5_) - (S 2 ))((S_S + ) (S 2 ))((S + S_ ) 2 - (5_5 + ) 2 ) A 
«S 2 ) 2 -(S_S + }<S + S-)) 2 

. 2eh ({S + S-} + (S 2 ))((S_S + ) + (SM((S+S-) - (S 2 )) 2 + ((S-S+) - (S, 2 )) 2 ] A 

((5 2 ) 2 -(5_5 + )(5 + 5_)) 2 H 

(17) 

Both the charge Hall conductivity and the spin Hall conductivity are inversely proportional 
to the absolute value of the chemical potential, and independent of the electron-impurity 
interaction strength and the concentration of scatterers. From Eqs. fTBf and (|1T[) we deduce 
that the charge and spin Hall effect should in principle be observable in graphene doped 
with magnetic impurities by changing the chemical potential close to the gap. As for the 
presence of SOI in graphene, we note that in recent experiments graphene is deposited on 
top of a Ni(lll) substrate. It is argued that the formation of a charge density gradient in the 
interface layer results in the experimentally observed large Rashba splitting in graphene [35]. 
Another aspect is that intercalated Au provides a « 100-fold enhancement of the spin-orbit 
splitting of graphene ir states [36]. In addition, as an open surface, graphene allows precise 
adatoms manipulations, as recently demonstrated experimentally [37]. Furthermore, when 
graphene is deposited on a ferromagnetic material such as Ni, the coupling between the 
Dirac fermion in graphene and the magnetic atom in the substrate arises inevitably due 
to roughness in the underlying substrate surface [38,39]. Motivated by these facts, we are 
hopeful that an possible experimental realization of charge and spin Hall effect in graphene 
in the presence of magnetic impurities, as uncovered here, should be possible. 

Summary. — We studied the charge and spin Hall effects in graphene with magnetic 
impurities. Using the Kubo formula, we obtained analytical expressions for the charge and 
spin Hall conductivities and concluded that both are generally finite in the presence of 
the SOI and magnetic impurities under zero external magnetic field. The charge Hall effect 
originates from a combined action of SOI and spin flip scattering at the magnetic impurities. 
The SOI in graphene results in transverse forces different for the two spin channels which 
yields a spin Hall current. On the other hand, the scattering from the magnetic impurities 
act as a spin-dependent barrier causing the imbalance of the charge accumulation at the 



p-5 



Ding et al. 



boundaries. The derived results for the charge and spin Hall effects should be observable by 
current technology [37] by changing the chemical potential close to the gap. 
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